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Abstract. The intermediate-frequency-peaked BL 
Lac object 3C 66A was detected in a flaring state 
by the Fermi-LXT and VERITAS observatories in 
October 2008. These data and follow-up observations 
at other wavelengths create a rich sample of light 
curves and a constraining spectral energy distribu- 
tion (SED). This is the first time that simultaneous 
observations at GeV and TeV energies were obtained 
for a flaring blazar. Results from these joint Fermi- 
LAT and VERITAS observations are presented in 
this paper. 
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I. Introduction 

Due to its significant optical and X-ray variability, 3C 
66A was classified as a BL Lac object by Maccagni et 
al. fjQ, an d given the location of its synchrotron peak 
(between 10 15 and 10 16 Hz), 3C 66A can be further 
sub-classified as an intermediate-frequency peaked BL 
Lac (IBL). BL Lacs are known for having very weak (if 
any) detectable emission lines, which makes the redshift 
determination quite difficult. The redshift of 3C 66A 
was reported as z — 0.44 by Miller et al. and 
independently by Lanzetta et al. [3|. However, both 
measurements rely on the measurement of a single line, 
and as pointed out by Bramel et al. [4|, the redshift of 
this object is still quite uncertain. 

Similarly to other blazars, the spectral energy distribu- 
tion (SED) of 3C 66A is known to have two pronounced 
peaks (see [5| for SED details), which suggests at least 
two different physical emission processes at play. The 
first peak (extending from radio to soft X-ray frequen- 
cies) is likely due to polarized synchrotron emission 
from high energy electrons, while different emission 
models (briefly introduced below) have been proposed 
to explain the second energy peak, which extends up to 
gamma-ray energies. It should be noted that the high- 
energy emission component dominates the bolometric 
luminosity of this source, a powerful demonstration of 
the extreme nature of this type of object. 

The models that have been proposed to explain 
gamma-ray emission in blazars can be roughly cat- 
egorized into leptonic and hadronic mechanisms, de- 
pending on whether the accelerated particles respon- 
sible for the gamma-ray emission are primarily elec- 
trons/positrons or protons. In leptonic models, high- 



energy electrons/positrons produce gamma rays via in- 
verse Compton scattering of low-energy photons. In 
synchrotron self-Compton (SSC) models, the same pop- 
ulation of electrons responsible for the observed gamma 
rays generates the low-energy photon field through syn- 
chrotron emission. In external Compton (EC) models the 
low-energy photons originate outside the jet. Possible 
sources of external photons include: accretion disk pho- 
tons radiated directly into the jet, accretion disk photons 
scattered by emission-line clouds or dust into the jet, or 
synchrotron radiation re-scattered back into the jet by 
broad-line emission clouds. In hadronic models gamma 
rays are produced by high-energy protons, either via 
proton synchrotron radiation, or via secondary (see [6] 
and references therein for a review of blazar gamma-ray 
emission processes). 

One of the main obstacles in the broadband study of 
gamma-ray blazars is the lack of simultaneity, or at least 
contemporaneity, of the data at the various wavelengths. 
At high energies the situation has been even more 
difficult due to the lack of objects that can be detected 
by MeV/GeV and TeV observatories in comparable time 
scales. Indeed, until recently most of our knowledge 
of blazars at gamma-ray energies was obtained from 
observations performed in two disjoint energy regimes: 
i) the high energy (HE) range (20 MeV< E < 10 
GeV), studied in the 1990s by EGRET Q, and ii) 
the very-high-energy (VHE) regime (E > 100 GeV) 
observed by ground-based instruments |8|. Remarkably, 
most EGRET blazars have not been detected by TeV 
telescopes; even those that are nearby and bright (9). 
Furthermore, blazars detected by EGRET at MeV/GeV 
energies are predominantly FSRQs, while blazars TeV 
blazars are BL Lacs. It is important to understand these 
observational differences since they are likely related 
to the physics of the AGN iflOl . or to the evolution of 
blazars over cosmic time fill . 

Fortunately, the new generation of gamma-ray instru- 
ments (AGILE, Fermi, HESS, MAGIC and VERITAS) 
is closing the gap between both energy regimes due 
to their improved sensitivities, leading us towards a 
deeper and more complete characterization of blazars as 
a high-energy source and as a population. An example of 
the successful synergy of space-borne and ground-based 
observatories was provided by the joint observations 
of 3C 66A by Fermi and VERITAS during its strong 
flare of October 2008. The strong flare was originally 
reported by VERITAS fl2l . and soon after, correlated 
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variability was detected in the Fermz-LAT data lfl3l . 
Follow-up observations were obtained at radio, optical 
and X-ray wavelengths in order to measure the flux and 
spectral variability of the source across the electromag- 
netic spectrum and to obtain a quasi-simultaneous SED. 
In this contribution we report preliminary results of this 
campaign, the resulting broadband spectrum and our first 
attempt to model this very constraining SED. 




Fig. 1 : Smoothed count map of the 3C 66A region as seen by Fermi 
between Sep 1 and Dec 31 2008 with E>100 MeV. The color bar has 
units of counts per pixel and the pixel dimensions are 0.1 x 0.1 degrees. 
The locations of 3C 66A and 3C 66B are shown as crosses and labeled 
"A" and "B" respectively. The location of PSR 0218+4232 is also 
indicated. The magenta circle represents the VERITAS localization of 
the VHE source (2h 22m 41.6s ±1.7s (stat) ±6.0s (sys) ; 43° 02' 35.5" 
±21" (stat) ±1'30" (sys)) as reported in Acciari et al. 1141 . The red 
circle represents the 95% error radius of the Fermi localization (RA, 
DEC) = (35.6541, 43.0402) and includes statistical error 0.042° and 
a conservative systematic error of 0.04° as discussed in Abdo et al. 
[20]. Although the Fermi localization clearly favors blazar 3C 66A as 
the gamma-ray source counterpart, some small contribution from radio 
galaxy 3C 66B 1 1 5 ] cannot be excluded. A more detailed analysis is 
ongoing to determine the likelihood of such a scenario. 

II. Observations and Data Analysis 

VERITAS: The Very Energetic Radiation Imaging 
Telescope Array System (VERITAS) is an array of 
four 12m diameter imaging atmospheric Cherenkov 
telescopes (IACTs) in southern Arizona, U.S.A. Ifl6l . 
VHE emission from 3C 66A was detected by VERITAS 
from 33 hours of observations between September and 
November 2008. A complete description of the VER- 
ITAS analysis and results is available in Acciari et al. 
fT4"l . It suffices to summarize here that a strong flare with 
day-scale variability was observed in October 2008 (as 
illustrated in the bottom panel of fig. 2) and that the 
observed spectrum is very steep (as shown in fig. 3), 
yielding an index of 4.1 ± 0.4 (stat) ± 0.6 (sys) when 
fitted to a power law. 

Fermi-hAT: The data from the Large Area Telescope 
(LAT; [17 1) have been analyzed by using ScienceTools 
v9rll. Only class-3 events (also called diffuse) were 



selected for this analysis because of their negligible 
background contamination and very good angular re- 
construction. In order to avoid Earth's albedo photons, a 
zenith angle <105 degrees cut in instrument coordinates 
was used. The diffuse emission from the Galaxy was 
modeled using a GALPROP model |18l] which has 
been refined with Fermi-LPR data taken during the first 
3 months of operation. The extragalactic diffuse and 
residual instrumental backgrounds have been modeled 
as an isotropic component and included in the fit (see 
l20l for a complete discussion of the treatment of 
these backgrounds). The data were analyzed with an 
unbinned maximum likelihood technique |fl9ll using the 
likelihood analysis software developed by the LAT team. 
Although 3C 66A was detected by EGRET as source 
3EG J0222+4253 ED, Kuiper et al. El showed with 
detailed spatial and timing analyses that this EGRET 
source actually consists of the superposition of 3C 66A 
and nearby millisecond pulsar PSR J02 18+4232. This 
interpretation of the EGRET data has been verified by 
ferm/-LAT, whose improved angular resolution permits 
the clear separation of the two sources as shown in 
Fig. 1. Furthermore, gamma-ray pulsations from PSR 
J02 18+4232 have been detected with high confidence 
by the Fermi team and will be reported elsewhere ||231 . 
More importantly for this analysis, the clear separation 
between the pulsar and the blazar allows us to account 
for each source independently in the maximum likeli- 
hood analysis and thus, permits an accurate determina- 
tion of the spectrum and localization of each source. 

Chandra: 3C66A was observed by the Chandra Ob- 
servatory on October 6th 2008 for a total of 37.6 ksec 
with the Advanced CCD Imaging Spectrometer (ACIS), 
covering the energy band between 0.3 and 10 keV. 
The source was observed in the continuous clocking 
(CC) mode to avoid pile-up effects. The standard data 
processed files provided by the Chandra X-ray Center 
(processed with software ver. 7.6.11.9 and calibration 
files (CALDB) ver. 3.5.0) have been used. 

Swift XRT and UVOT: Following the VERITAS 
detection of VHE emission from 3C 66A, Target of 
Opportunity (ToO) observations of 3C66A with Swift 
were obtained for a total duration of ~ lOksec. Swift 
comprises an UV-Optical telescope UVOT, an X-ray 
telescope (XRT) and a Burst Alert Telescope |24l . 
Data reduction and calibration of the XRT data were 
performed with HEASoft v6.5 standard tools. All XRT 
data presented here were taken in Photon Counting (PC) 
model with negligible pile-up effects. The X-ray spectra 
of each observation were fit with an absorbed power 
law using a fixed Galactic column density from Dickey 
& Lockman |25l . which gave good \ 2 values for all 
observations. UVOT observations were taken over six 
bands (V, B, U, Wl, M2 and W2) that were calibrated 
using standard techniques l26l . Il27ll and redshift cor- 
rected using a nominal value of z = 0.44. 

MDM Optical Observations: 3C 66A was observed 
at the 1.3m telescope of the MDM Observatory during 
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Fig. 2: 3C 66A light curves from (top to bottom): PAIRITEL, MDM, 
Swift XRT, Fermi and VERITAS. The Fermi panel includes light 
curves for E > 200 MeV and E > 1 GeV using time bins with width 
equal to 3 days. It should be noted that the flare observed by VERITAS 
coincides with the period of time when Fermi detected the source at 
high energies (> 1 GeV). 



the nights of Oct. 6- 10, 2008 UT. A total of 290 science 
frames in U, B, V, R, and I bands (58 each) were taken 
throughout the entire visibility period (approx. 4:30 - 
10:00 UT) during each night. 

PAIRITEL Observations: Near-infrared observa- 
tions in the J, H and K s were obtained with the 1.3m 
Peters Automated Infrared Imaging Telescope located 
at the Fred Lawrence Whipple Observatory. As a fully 
robotic observatory, PAIRITEL relies on a software 
system that manages the telescope and reduces data as it 
is gathered. See Bloom et al. [28 1 for a full description 
of the data reduction procedure. 

Additional multiwavelength observations at radio, op- 
tical and X-ray wavelengths will be included in a future 
paper. 

III. Discussion 

The measured SED of 3C66A in October 2008 is 
presented in Fig. 3. We have modeled the SED using the 
code of Bottcher & Chiang J29]. In this model, a power- 
law distribution of ultrarelativistic electrons and/or pairs 
with lower and upper energy cutoffs at 7 m ; n and j max , 
respectively, and power-law index q is injected into a 
spherical region of co-moving radius Rb- The injection 
rate is normalized to an injection luminosity L e , which 
is a free input parameter of the model. A temporary 
equilibrium between particle injection, radiative cooling 
due to synchrotron and Compton losses, and particle 
escape on a time scale t csc = r] csc Rb /c is assumed. 
Both the internal synchrotron photon field (SSC) and 
external photon sources (EC) are considered as targets 
for Compton scattering. The emission region is moving 
with a bulk Lorentz factor V along the jet. To reduce 
the number of free parameters, we assume that the 
jet is oriented with respect to the line of sight at the 
superluminal angle so that the Doppler factor is equal 
to D = (r [1 — f3 cos^obs]) -1 where o bs lS the angle 
of the jet with respect to the line of sight. We use the 
extragalactic -background-light model of Franceschini et 
al. ||30 l and the adopted value of z = 0.444 for our 
model, but keep in mind the caveat that the results 
might be significantly different if the source is at a very 
different redshift. 



Parameter 


Symbol 


Value 


e~ Injection Luminosity 


L e 


1.26 X 10 4b erg s" 1 


Low-energy cut-off 


7min 


5 X 10 3 


High-energy cut-off 


7max 


8 x 10 4 


e~ injection index 


q 


2.5 


Plasmoid radius 


Rb 


2 X 10 16 cm 


Bulk Lorentz factor V 


r 


30 


Observing Angle 


#obs 


1.9 deg 


Co-moving magnetic field 


B 


0.3 G 


e~ escape time scale 




100 


Ext. radiation energy density 


Mcxt 


6.6 X 10" 6 erg cm" 3 


Ext. radiation peak frequency 


^cxt 


1.5 x 10 14 Hz 



TABLE I: EC+SSC model fit parameters 

Most TeV blazars known to date are high-frequency 
peaked BL Lacs (HBLs), whose SEDs can often be fit 
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Fig. 3: SED of 3C 66A including Fermi-LAT and VERITAS 
data. The Fermi spectrum was calculated for the short time interval 
Oct 4 - 6 (MID 54743-54745) where most of the multiwavelength 
observations took place and for a longer time interval corresponding 
to the VERITAS dark period where VHE emission from 3C 66A was 
detected (MID 54734 - 54749). The Fermi flux from Oct 4 - 6 is 
higher F(E>100MeV) = 4.9 ± 0.9 x 10~ 7 cm- 2 s _1 than the dark 
period flux of 3.6 ± 0.4 x 10 _7 cm _2 s — 1 . In both cases the Fermi 
spectrum is quite hard and can be described by a power law with index 
1.84 ± 0.11 and 1.87 ± 0.06 in the Oct 4-6 and dark period intervals 
respectively, thus, no spectral evolution is evident in the Fermi data. 
Furthermore, it should be noted that there are not enough photons at 
high-energies in the Fermi data (6 photons with E > 10 GeV for the 
Oct 4-6 period) to confirm the spectral steepening around 10 GeV 
indicated by the model. 

satisfactorily with pure SSC models. Therefore, we first 
attempted a pure SSC model to fit the near-IR through 
VHE gamma-ray SED. However, we found that such a 
model (dashed line in Fig. 3) grossly overproduces the 
gamma-ray flux at GeV energies and requires a magnetic 
field far below the value expected from equipartition. 
A better fit to the entire SED could be achieved by 
including an additional external radiation field, modeled 
as a thermal radiation field with a peak frequency at 
^ext = 1-5 x 10 14 Hz. and a radiation density of u cx t = 
6.6 x 1(T 6 erg cm" 3 . Nevertheless, the EC+SSC model 
(solid line in Fig. 3) still requires a sub-equipartition 
magnetic field (ub/u b ~ 0.1) (the full set of parameters 
used for the fit shown in Fig. 3 is listed in table 1). 
These preliminary results agree with the expectation 
that for IBLs the higher ambient luminosity at optical 
wavelengths plays a fundamental role in providing a 
seed population of soft photons for inverse Compton 
scattering. 

IV. Conclusion and Future Plans 

Multiwavelength observations of 3C 66A were carried 
out following the gamma-ray outburst detected by the 
VERITAS and Fermi observatories in October 2008. 
This marks the first occasion that a gamma-ray flare is 
detected by GeV and TeV instruments in comparable 
time scales. The SED constructed with these and other 



contemporaneous multiwavelength observations consti- 
tutes a fertile ground for the modeling of high-energy 
emission in gamma-ray blazars. A joint paper to further 
address this and other topics is in preparation by the 
Fermi-LAT and VERITAS collaborations. 
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